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Abstract. We present a new mechanism to explain the 
frequently observed and thus certainly permanent warping 
of spiral galaxies. We consider the possibility of non-linear 
coupling between the spiral wave of the galaxy and two 
warp waves, such that the former, which is linearly unsta- 
ble and extracts energy and angular momentum from the 
inner regions of the galactic disk, can continuously feed the 
latter. We derive an expression for the coupling coefficient 
in the WKB approximation. We show that the coupling is 
too weak in the stellar disk, except at the Outer Lindblad 
Resonance where the spiral slows down and is efficiently 
coupled to warp waves. There, the spiral can be almost 
completely converted into "transmitted" warps, which we 
can observe in HI, and a "reflected" one, which we can 
observe as a corrugation. Our mechanism reproduces the 
observed amplitudes of the warp and of the corrugation, 
and might explain related phenomena such as the behav- 
ior of the line of nodes of the warp. Furthermore we show 
that the energy and momentum fluxes of observed spirals 
and warps are of the same order of magnitude, adding a 
strong point in favor of this model. 

Key words: Galaxies : warping, spiral waves, Lindblad 
resonances - Nonlinear coupling 



1. Introduction 

Spiral galaxies often show a warping of their external 
parts, observed in HI. Such warps appear like an "S" or an 
integral sign for galaxies observed edge-on (Sancisi, 1976). 
This phenomenon, known since 1957, has led to theo- 
retical difficulties. Warps are bending waves, and we cor- 
rectly understand their propagation (see Hunter, 1969a 
for a dispersion relation in an infinitely thin disk, and see 
Nelson 1976a and 1976b, Papaloizou and Lin 1995, Mas- 
set and Tagger 1995 for the dispersion relation taking into 
account finite thickness and compressional effects), but we 
do not know the mechanism responsible for their excita- 
tion. It cannot be systematically justified by tidal effects 
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(Hunter and Toomre 1969), since we observe warping of 
very isolated spiral galaxies. It cannot be explained by a 
temporary excitation, since a warp propagates radially in 
the disk is not reflected at its edge, so that it should disap- 
pear over a few galactic years (Hunter, 1969b). Attempts 
have been made to connect the existence of warps and the 
properties of haloes (Sparke, 1984a, Sparke and Casertano 
1988, Hofner and Sparke 1994). The mechanism proposed 
by Sparke and Casertano (1988) implies an ad hoc mis- 
alignment angle between the normal axis of the galactic 
plane and the "polar" axis of the halo. This mechanism 
has recently encountered self-consistency difficulties (Du- 
binski and Kuijken 1995). Furthermore, it relics on the 
dubious hypothesis of a rigid and unresponsive halo. Sev- 
eral explanations have also been proposed, such as the in- 
fall of primordial matter (Binney, 1992), the effect of the 
dynamical pressure of the intergalactic medium, etc. . . . 
For a review of all the proposed mechanisms, see Binney 
(1992). 

Binney (1978, 1981) has considered the possibility of 
a resonant coupling between the vertical motion of a star 
and the variation of the galactocentric force, due to a halo 
or a bar. He concluded that a bar could be responsible for 
the observed warps and corrugations. Sparke (1984b) has 
also explored this possibility, and considered the growth 
of a warp from a bar or a triaxial halo. She found that 
a bar was unlikely to be responsible for a warp, but she 
emphasized that a triaxial halo could quite well reproduce 
observed warps. 

We consider here another mechanism, the non-linear 
coupling between a spiral wave and two warp waves. Non- 
linear coupling between spirals and bars has already been 
found (Tagger et al. 1987 and Sygnet et al. 1988) to pro- 
vide a convincing explanation for certain behaviors ob- 
served in numerical simulations (ScUwood 1985), or ob- 
served in Fourier Transforms of pictures of face-on galax- 
ies. The relative amplitude of spiral or warp waves (the 
ratio of the perturbed potentials to the axisymmetric one) 
is about '--^ 0.1 to ~ 0.3 (see Strom et al. 1976). Non-linear 
coupling involves terms of second order in the perturbed 
potential, while linear propagation is described by first 
order terms. Classically one would thus believe that non- 
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linear coupling is weak at such small relative amplitudes. 
However in the above-mentioned works it was found that 
the presence of resonances could make the coupling much 
more efficient if the wave frequencies are such that their 
resonances (i.e. the corotation of one wave and a Lind- 
blad resonance of another one) coincide. At this radius the 
non-linear terms become comparable with the linear ones, 
so that the waves can very efficiently exchange energy and 
angular momentum. Indeed in Sellwood's (1985) simula- 
tions, as discussed by Tagger et al. (1987) and Sygnet et al. 
(1988), an "inner" spiral or bar wave, as it reaches its coro- 
tation radius, transfers the energy and angular momentum 
extracted from the inner parts of the disk to an "outer" 
one whose ILR lies at the same radius, and which will 
transfer them further out, and ultimately deposit them at 
its OLR. In this process the energy and momentum are 
thus transferred much farther radially than they would 
have been by a single wave, limited in its radial extent by 
the peaked rotation profile. 

We will show here that a similar mechanism, now in- 
volving one spiral and two warp waves is not only possible 
(by the "selection rules" associated with their parity and 
wavenumbers), but also very efficient if the same coinci- 
dence of resonances occurs. This allows the spiral wave, 
as it reaches its OLR (and from linear theory deposits the 
energy and momentum extracted from the inner regions 
of the disk) to transfer them to the warps which will carry 
them further out. 

Unhkc Tagger et al. (1987) and Sygnet et al. (1988), we 
will throughout this paper restrict our analysis to gaseous 
rather than stellar disks, described from hydrodynamics 
rather than from the Vlasov equation. The reason is that 
our interest here lies mainly in the excitation of the warps, 
which propagate essentially in the gas (indeed the outer 
warp is observed in HI, and the corrugation is most likely 
(Florido et al. 1991) due to the motion of the gaseous 
component of the galactic disk). On the other hand, the 
spiral wave propagates in the stellar as well as the gaseous 
disk. The difference is important only in the immediate 
vicinity of Lindblad resonances, where the spiral wave is 
absorbed; as a consequence, its group velocity vanishes 
at the resonances. Since the group velocity of the waves 
will appear as an important parameter, we will choose to 
keep the analytic coupling coefficient derived from the hy- 
drodynamic analysis, but we will introduce, for the spiral 
density wave, the group velocity of a stellar spiral. From 
the physics involved this will appear as a reasonable ap- 
proximation; furthermore it should only underestimate the 
coupling efficiency, since it does not include the resonant 
stellar motions near the resonance. 

On the other hand, we will show that non-linear cou- 
pling is efficient only in a narrow annulus close to the OLR 
of the spiral, over a scale length similar to the one of Lan- 
dau damping. We will thus conclude that the two processes 
are in direct competition, with the spiral transferring its 
energy and momentum, in part to the stars by Landau 



damping, and in part to the warps which will transfer 
them further outward, the exact repartition between these 
mechanisms presumably depending on detailed character- 
istics of the galactic disk. 

The paper is organized as follows: in a first part we 
will introduce the notations, and the selection rules rela- 
tive to the coupling. In a second part, we will derive the 
coupling coefficient from the hydrodynamic equations ex- 
panded to second order in the perturbed quantities, and 
we will try and simplify it. In a third part, we will analyze 
the efficiency of the coupling, together with the locations 
where it may occur. In the last sections we will compare 
our predictions to the observations, and we will propose 
some possible observational tests of our mechanism. 

Some of the computations are tedious and lengthy. For 
the sake of clarity, they are developed in appendices, so as 
to retain in the main text only the principal results and 
the physical discussions. 

2. Formalism and selection rules 

2.1. Notations 

We develop all our computations in the well-known shear- 
ing sheet approximation, which consists in rectifying a nar- 
row annulus around the corotation radius of the spiral 
wave into a Cartesian slab. The results given by an exact 
computation taking into account the cylindrical geome- 
try of the galaxy would differ by some metric coefhcients, 
but would not differ physically from the shearing sheet 
predictions, so that the main conclusions would remain 
valid. We call x the radial coordinate, with its origin at 
the corotation of the spiral, and oriented outward. We 
call y the azimuthal coordinate, oriented in the direction 
of the rotation, and z the vertical one, so as to construct 
a right-oriented frame xyz. 

The hydrodynamic quantities are the density p, the 
speed components {U, V, W), and the gravitational poten- 
tial (p. We assume that the disk is isothermal, with a uni- 
form temperature, in order to avoid unnecessary complex- 
ity. Thus we can write: P ~ a^p, where P is the pressure 
and a the isothermal sound speed. 

The perturbed quantities are denoted with subscripts 
1, 2 or S" referring to the wave involved (one of the two 
warps or the spiral wave), and equilibrium quantities with 
a subscript 0. 

The epicyclic frequency is denoted by k, the rotation 
frequency by fi, Oort's first constant by ^ = O'r/2, and 
we have the relation: 

K^ = An{n + A) 

We call /i the vertical characteristic frequency of the 
disk. We have: 



M 



2n^ 



(see Hunter and Toomre, 1969). This frequency plays 
the same role for bending waves as k does for spiral waves; 
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it should not be confused with the frequency v^ of the ver- 
tical motion of individual particles, which is usually much 
higher. Indeed individual particles move vertically in the 
potential well of the disk, giving them the high frequency 
Vz ■ On the other hand, when one considers motions of the 
whole disk, the potential well moves together with the disk 
and exerts no restoring force. This leaves only a weaker 
restoring force, giving the vertical frequency [i as discussed 
e.g. in Hunter and Toomre (1969) or Masset and Tagger 
(1995). 

For each of the three waves we consider, we denote by w 
its frequency in the galactocentric frame, m its azimuthal 
wavenumber (jn = 2 for a two-armed spiral, and m = 1 
for an "integral-sign" warp), w = w — m^(r) its frequency 
in a local frame rotating with the matter. In the peculiar 
case of the shearing sheet, we note ky = mjr where r is 
the distance to the galactic center, and ky the azimuthal 
wavenumber. 

We will have to perform part of our computations in 
the WKB approximation, i.e. assume that the radial 
wavevector varies relatively weakly over one radial wave- 
length; in practice, in the shearing sheet, this reduces to 
the classical "tightly wound" approximation, ky <^ kx. We 

will also use q = Jk^ + ky , the modulus of the "horizon- 
tal" wavevector. In the tightly wound approximation one 
has q ~ \kx\. 

We also introduce s — p/ Pq, and various integrated 
quantities: 



+ 00 



p{z)dz 



the perturbed surface density, 
S = / po{z)dz 



the equilibrium surface density, and: 



Z = - 



+ 00 



zp{z)dz 



the mean vertical deviation of a column of matter from 
the midplane, under the influence of a warp. 

We call H the characteristic thickness of the disk. The 
vertical density profile in the disk is taken to be consistent 
(see Masset and Tagger, 1995), i.e. it must fulfill simulta- 
neously the Poisson equation and the hydrostatic equi- 
librium equation, with the additional condition that the 
radial derivatives of the equilibrium potential do not de- 
pend on z throughout the disk thickness: this condition 
results from the hypothesis of a disk which is geometri- 
cally thin, if <C r, although we do resolve vertically the 
perturbed quantities along the vertical direction. 

We summarize most of our notations in figures |l| and 
g and in table |^. 
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h ig. 1. This figure summarizes our main notations relative to per- 
turbed velocity and disk thickness. 
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X* ig. 2. This figure summarizes our notations relative to waves (spiral 
or warps). The column of matter has a section unity. At equilibrium 
(or more generally without spiral wave) it contains a mass S, and when 
perturbed by a spiral wave it contains the mass S + cr. 



Symbol 


corresponding quantity 


a 


Sound speed 


Q. 


Equilibrium angular velocity of matter 


A 


Oort's first constant. A — l/2rdQ,/dr 


K 


Epicyclic frequency, k. — 4S1(0 + A) 


P 


"Global vertical frequency", fi — 20 — k 


LO 


Wave frequency in galactocentric frame. 


Q 


Wave frequency in rotating frame. 





Perturbed potential. 


p 


Perturbed density. 


s 


Relative perturbed density (s — p/po). 


Q 


Wave- vector modulus (g — {k^ -\- A;^)"'"'^) 



J.a.ble 1. In this table we summarize the main frequencies of our 
system, and some other quantities related to warps or spirals. 
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2.2. Notion of coupling and selection rules 

2.2.1. Mode coupling 

As discussed in the introduction, we consider the coupling 
between a spiral wave and two warp waves. This coupling 
is necessarily a non-linear mechanism. In a linear analy- 
sis, each wave can be studied independently, by projecting 
perturbed quantities onto the frequency and the wavevec- 
tor of each wave. In particular, combining the hydrody- 
namic (continuity and Euler) equations and the Poisson 
equation leads to the dispersion relation which reads, in 
the WKB limit: 



a)| = K^ - 27rGE(7s + a^ql 




for a spiral and: 




~ 2 
,-.2 2 1 o__/T^„,, 1 ^W 


— fl2„2 


^W ^ h^ ' /jTTI^ZjqiY + ~2 


2« 9W 



for a warp (see Masset and Tagger 1995, Papaloizou and 
Lin 1995). 

Here we expand the hydrodynamic equations to sec- 
ond order in the perturbed quantities. Then the waves 
do not evolve independently anymore and, provided that 
they satisfy selection rules which will be discussed below, 
they can exchange energy and angular momentum. 

Thus the problem we address can be described as fol- 
lows: 

We consider a spiral wave, with a given flux (its en- 
ergy density multiplied by its group velocity), traveling 
outward, from its corotation to the Holmberg radius of 
the galaxy. We will not discuss here how this spiral has 
been excited and wc do not try to justify its amplitude, 
but rather take it as an observational fact. 

As it travels radially the spiral interacts with warps 
which are always present, at all frequencies, at a noise 
level associated with supernovae explosions, remote tidal 
excitations, etc. 

Our goal is the following: can wc find conditions such 
that, although the warps are initially at this low noise 
level, they can be non-linearly coupled to the spiral so ef- 
ficiently that they absorb a sizable fraction of its energy 
and momentum fiux ? And can these conditions be met 
commonly enough to explain the frequent (one might even 
say general) occurrence of warps and their main observa- 
tional properties ? 

We will first discuss the conditions, known as "selec- 
tion rules" , for the warps to be coupled to the spiral. 

2.2.2. Selection Rules 

The linearized set of equations governing the wave behav- 
ior (continuity, Euler and Poisson) is homogeneous and 
even in z, so that any field of perturbations to the equi- 
librium state of the disk can be considered as composed 
of two independent parts: 



— Perturbations whose perturbed density is even in z, 
and thus, due to hydrodynamic equations, whose per- 
turbed quantities are all even in z except W, which is 
odd. These perturbations are spiral waves, since they 
imply a perturbed density (a non- vanishing a), and a 
vanishing Z (so they do not raise the mid-plane of the 
disk). 

— Perturbations whose perturbed density is odd in z, and 
thus whose perturbed quantities are all odd in z, ex- 
cept W, which is even. These perturbations are warps, 
since they imply a vanishing integrated perturbed den- 
sity and a non- vanishing Z, which means that they 
involve a global motion of the mid-plane of the disk. 

Furthermore, since the coefficients of the equations do 
not depend on time or azimuthal angle, Fourier analysis 
allows to separate solutions identified by their frequency 
and azimuthal wavenumber. On the other hand, since the 
coefficient do depend on the radius, Fourier analysis in r 
(or X in the shearing sheet) does not allow the definition of 
a radial wavenumber except in the WKB approximation. 
We will return to this below. 

So let us consider a spiral wave for which we write: 

^s oc e ^ -^ ^ ' 

and two warp waves with: 

Ci ex e'("i''-'^i*) 

and 

^2 oc e'(™2^-'^2t) 

where ^ represents any perturbed quantity. Hydro- 
dynamic equations written to second order in perturbed 
quantities will contain terms involving the products ^1^2, 
ei?2*, er6, etc. ^ 

Let us consider for instance the term ^1^2 (rules for 
the other products are derived in a similar manner). Its 
behavior with time and azimuthal angle is: 

c^c^ ex e*[(™i+™2)s-('-'i+'^2)t] 

More physically this product, which appears from such 
terms as yi.VV2 or piV02, can be interpreted as a beat 
wave. Fourier analysis in t and 9 will give a contribution 
from this term, at frequency loi + loi and wavenumber 
m\ -\- TO2; thus it can interact with the spiral wave if: 

0^1+6^2= ^S 

and 

m\ -}- TO2 — rris 

As discussed above, the equations for the spiral wave 
are derived from the part of the perturbed quantities 
which is even in z. Our third selection rule is thus that 
the product ^1^2 be even in z. 



F. Massct & M. Tagger: Excitation of warps by spiral waves in galaxies by non-linear coupling 



If we were in a radially homogeneous system, so that 
waves could also be separated by their radial wavenumber, 
a fourth selection rule would be: 

kxi + kx2 = kxs- 

Since this is not the case, we will write a coupling coef- 
ficient which depends on x. This coefficient would vanish 
by radial Fourier analysis if the waves had well-defined ra- 
dial wavenumbers, unless these wavenumbers obeyed this 
fourth selection rule. We will rather find here that the co- 
efficient varies rapidly with x and gives a coupling very 
localized in a narrow radial region. This is not a surprise 
since it also occurred in the interpretation by Tagger et 
al. (1987) and Sygnet et al. (1988), in terms of non-Hnear 
coupling between spiral waves, of the numerical results of 
Sellwood (1985). The localization of the coupling, as will 
be described later, is associated with the presence of the 
Lindblad resonances. We will find that, in practice, it re- 
sults in an impulsive-like generation of the warp waves, in 
a sense that will be discussed in section 4.2.2. 

Let us mention here that we have undertaken Particle- 
Mesh numerical simulations in order to confirm this analy- 
sis. Preliminary results in 2D, with initial conditions simi- 
lar to Sellwood (1985), do show the to = and m = 4 
spiral waves, at the exact frequencies and radial loca- 
tion predicted by Tagger et al. (1987) and Sygnet et al. 
(1988). The same absence of a selection rule for the radial 
wavenumber is observed. These results will be reported 
elsewhere, and in a second step the simulations will be 
applied in 3D to the physics described in this paper. 

A radial selection rule would allow us to get a very 
simple expression for the efficiency of non-linear coupling 
-i.e. for wave triplets that do not obey it, and 1 for 
triplets that do. Here we will have to rely on a more deli- 
cate integration of the coupling term over the radial extent 
where it acts. This will be done in section 4.2. In particu- 
lar in 4.2.2, we will show that in the vicinity of the OLR of 
the spiral, even though its WKB radial wavenumber is di- 
vergent, the main localization of the coupling comes from 
an (integrable) divergence of the coupling coefficient. 

Let us summarize the selection rules: 

~ We first have the condition to-i -I- Tn^ — rns. This is 
obviously fulfilled by the most frequent warps (which 
have TO = 1) and spirals (with to = 2). One should 
note here that, since we are computing with complex 
numbers but dealing with real quantities, each pertur- 
bation is associated with its complex conjugate, with 
wavenumber —m and frequency —uo*. Thus we find a 
contribution to the warp 1 by the coupling of the spi- 
ral with the complex conjugate of warp 2, i.e. from 
products of the form S,s£,2i 6tc. 

— We then have the parity condition, which is obviously 
fulfilled by odd warps and an even spiral. 

— Finally the frequency selection rule: 

UJl+ UJ2— ^S 



gives us in principle an infinite choice of pairs of warp 
waves, since the latter can be presumed to form a con- 
tinuous spectrum. One of our main tasks in the present 
work will be to determine which pair is preferentially 
coupled to the spiral. Our assumption of coupling be- 
tween only three waves will be found valid when we 
find that actually one such pair is strongly favored, so 
that all the others can be neglected. 
The frequency condition can also be written as: 

a)i + a;2 = ti^i — rj + 0^2 — f^ = tJs — 2f7 = ujs 

showing that it is also true in the rotating frame at 
any radius, if it is true anywhere. 
A second remark concerning this selection rule is that 
Lo can have an imaginary part. One is thus restricted 
to two possible choices : either working in a two-time- 
scales approximation, where the fast scale is that of the 
oscillations ((-^^^), while the slow one is both that of 
linear growth or damping and of non- linear evolution; 
or simplify the problem by only looking for permanent 
regimes, keeping uj real. We will retain the second pos- 
sibility, and in fact consider the non-linear evolution of 
the waves as a function of x as they travel radially, as- 
suming that the inner part of the disk feeds the region 
we consider with spiral waves at a constant amplitude. 
This is sufficient for our main goal, which is to show 
that non- linear coupling is indeed possible and efficient 
to generate warps. On the other hand one should keep 
in mind that the permanent regime we will find may 
very well be unstable, since it is well known that mode 
coupling (in the classical case of a homogeneous sys- 
tem, much simpler than the one we consider) can lead 
to any type of complex time behavior, e.g. limit cycles 
or even strange attractors. 

3. The coupling coefficient 

Once we have the set of selection rules, we can derive an 
expression of the coupling coefficient between the spiral 
and the two warps. In a linear analysis, each wave propa- 
gates independently and is thus subject to a conservation 
law of the form: 

dtE + -dr{rc„E) = 
r 

where 9„ stands for the partial derivative with respect to 

any quantity u, E is, the energy density of the wave, Cg its 

group velocity. This relation simply means that globally 

the energy of the wave is conserved and advected at the 

group velocity Cg — duj/dkx (Mark, 1974). The r and 1/r 

factors in the spatial derivative come from the cylindrical 

geometry of the problem. In the framework of the shearing 

sheet we neglect them and we may write: 

dtE = dtE + dxicgE) = 

The vanishing right-hand side results from the fact 
that the wave does not exchange energy with other waves 
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or with the particles, and thus is neither amplified nor ab- 
sorbed. Mode coupling introduces in the right-hand side 
a new term describing the energy exchange between the 
waves. Since E is quadratic in the perturbed amplitudes, 
and mode coupling corresponds to going one order fur- 
ther in an expansion in the perturbed amplitudes, this 
new source term will be of third order. Its derivation is 
lengthy and technical, and we give it in separate appen- 
dices. We summarize it in the next sections. 

3.1. First step : exact expression of the coupling 

The first step consists in deriving an exact expression for 
this coupling coefficient. To do this we integrate the hy- 
drodynamic equations over z after transforming their lin- 
earized parts to write them as variational forms. After a 
few transformations, we find that the total time derivative 
of a variational form, which involves the perturbed quanti- 
ties associated with each wave (and which we will interpret 
as its energy density) , is equal to a sum of terms involving 
the three waves, which we will interpret as the coupling 
part, i.e. the energy exchanged between the waves. This 
derivation is done in appendix |A| for the case of one of the 
warp waves. Let us just mention that this first step has 
the following features: 

— The only approximations we make are that the disk 
is isothermal with a uniform temperature, and that 
the horizontal velocities in each wave, U and V , are in 
quadrature, a reasonable assumption which is asymp- 
totically true in the limit of the WKB regime. 

— We fully resolve the thickness of the disk. Thus, the 
coupling coefficient involves integrals over z of prod- 
ucts of perturbed quantities related to each wave. 

— For the sake of compactness, the derivation is made in 
tensor formalism. This allows us to have a restricted 
number a terms (in fact eight terms in the coupling 
coefhcient), although at this step of the derivation we 
do not make any hypothesis on the perturbed motions 
associated with a spiral or a warp. 

Finally, at the end of this primary step, we obtain for 
the time evolution of the energy density of warp 1: 



dt 



(ViVi*+ assist) +Qst0i+c 
Jpoa^V,\ssd's*^^ 



(1) 



a'sld^,{psV^ 



Poa^Vrs^d'ss 



a'.s*A{p;V,^) 



where d/dt means d/dt + dx{cg. ). The L.H.S. 

appears as the total derivative of the energy density of 
warp 1, as expected, and the R.H.S. represents the cou- 
pling term, since each of the integrals involve the per- 
turbed quantities of the other waves (spiral and warp 2). 

3.2. Second step : analytic result in the WKB 
approximation 

Equation (|l]) is too complex to be used directly. We show 
in appendix O how it can be simplified by expanding the 
implicit sums, and explicitly writing the perturbed quan- 
tities associated with the warps or the spiral, and then by 
using in the evaluation of this term the eigenvectors (i.e. 
the various components of the perturbation) derived from 
the linear analysis. We make some assumptions: 

— We assume that the spiral involves no vertical motions, 
i.e. Ws = 0, and that horizontal motions do not de- 
pend on z, i.e. dzUs = dzVs — 0. These two assump- 
tions lead to an eigenvector which is the one of the 
infinitely thin disk approximation, although we verti- 
cally resolve the disk. 

— Symmetrically, we assume that the vertical speed in a 
warp is independent of z, i.e. : d^Wi = dzW2 = 0, 
so that the perturbed density is obtained by a simple 
vertical translation of the equilibrium density profile. 
This is consistent only if there is no horizontal motions. 
On the other hand, the compressibility of the gas disk 
does introduce horizontal motions in the warp wave 
(Masset and Tagger, 1995); thus, in order not to loose 
a possible coupling through these horizontal motions, 
we retain in the expression of the coupling coefficient 
the horizontal velocities due to the warps. This may be 
important near the Lindblad resonances of the warp, 
where horizontal motions can become dominant. 

At the end of this second step, we get a new expression for 
dEi/dt (where Ei is the energy density of warp 1) where 
we have performed the integrals over z, and which can be 
written: 



dt 



El — \\Zi\Zi2(y 



(2) 



+ 1 ^id^ipsvi*) + 1 4>\^^[p*1yh) + c.c. 



where Ai is a long expression which docs not need to 
be reproduced here, and which depends on x. Equation (|2|) 
shows that warp 1 is coupled to the product of the ampli- 
tude of both warps (Zi and Z-i) and to the amplitude of 
the spiral (cr). In particular, we see that if the warp does 
not exist at all in the beginning (i.e. it has a vanishing 
amplitude), it will never grow. Thus it has to pre-exist 
at some noise level in order to be allowed to couple with 
another warp and the spiral. We see here an important 
difference with the harmonic generation by a single wave, 
since in that case the harmonic will be generated even it 
it does not pre-exist. Here we are rather concerned with 
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the production of warps as "sub-harmonics" of the spi- 
ral. They cannot be spontaneously created, as this would 
correspond to a breaking of azimuthal symmetry. 

Mathematically, we see that the L.H.S. of the expres- 
sion above is second order with respect to the perturba- 
tions amplitudes, and that the R.H.S. is third order. In 
fact, the L.H.S. comes from the linear part of hydrody- 
namic equations, and has been transformed into a sec- 
ond order expression (a variational form), and the R.H.S. , 
which comes from the non-linear part of the equations, 
has simultaneously been converted into a third order ex- 
pression. 

In the same manner we can derive, by swapping indices 
1 and 2, an expression for the time evolution of the energy 
density of warp 2: 



al. (1987) and Sygnet et al. (1988). In that case the cou- 
pling coefficients, obtained from kinetic theory, were found 
to involve two resonant denominators, corresponding to 
the resonances of two of the waves; this made the cou- 
pling very efficient if the two denominators vanished at 
the same radius. We do not find such denominators here, 
because the warp vertical velocity (represented by Z in 
equation 0) does not diverge at the Lindblad resonance, 
while the horizontal velocities associated with the spiral 
do. On the other hand we will recover here a similar prop- 
erty of strong coupling close to the resonances, because 
the group velocity of the waves becomes small (it goes to 
zero in the asymptotic limit), so that the waves can non- 
linearly interact for a long time. This will be discussed in 
more details below. 



—r-E2 — \2ZiZ2(J 

at 

In order to close our set of equations, we also have 
to follow the behavior of the spiral. In the absence of any 
coupling to other waves, for reasons of global conservation 
of energy, the time evolution equation of the spiral is given 
by: 



dt 



Es = -(Al + \2)ZiZ2(y 



We choose to neglect coupling to other waves here since 
they either belong to the dynamics of the spiral itself {e.g. 
generation of m = 4 harmonics), and are irrelevant here 
since we take the spiral as an observational fact, or involve 
other warps: but, as mentioned above, we will find below 
that one pair of warps is preferentially driven by the spiral. 

3.3. Final step : the simplification 

The third and final step is to give an estimate of the cou- 
pling coefficients A. This is done in appendix O. Simplify- 
ing this coefficient implies much discussion on the physics 
involved, e.g. the behavior of resonant terms near the 
Lindblad resonances, or the order of magnitude of ratios 
of characteristic frequencies for a realistic galactic disk, 
etc. These discussions are fully developed in appendix ^ 
and lead us to the result that the coupling coefficient, any- 
where in the disk, is always of the order of vIGjs, i.e. : 



dt 



El ^ V^UJsZlZ2(7 



(3) 



where Vz represents the frequency of vertical oscilla- 
tions of a test particle in the rest potential of the disk 
(it must not be confused with fi, which is the frequency 
of global oscillations of the galactic plane, achieved when 
the whole disk - stars and gas- is coherently moved up 
and down). We see that this coefficient does not contain 
any dependency on 1 and 2, so that Ai ~ A2 '^ ^1'^- 

An important difference occurs here with the case of 
coupling between spirals or bars, analyzed by Tagger et 



4. The coupling efficiency 

We have to solve, or at least analyze, the behavior of a sys- 
tem of coupled differential equations which is non-linear 
(since the right-hand-sides are products of the unknowns). 
Such a system can have subtle and varied solutions, such 
as limit cycles or strange attractors (chaotic behavior). 
This goes far beyond the scope of this paper and we will 
restrict ourselves to looking for stationary solutions, which 
might prove to be unstable but still will inform us on the 
efficiency of the coupling mechanism. Thus we suppress 
the partial t derivatives in the total derivative d/di, mak- 
ing solution much simpler. 

J^.l. Search for a stationary solution 

We can write the system of three coupled equations in a 
slightly different manner, by introducing the energy fluxes 
(Hi = Cg.Ei) of the waves, and by relating density energies 
and amplitudes through constants Ki derived in appendix 

u 

E, = T.KfZf for i = 1,2 

Es^Kl^ 

(see for the derivation of Ki); the energy densities are 
proportional to the square of the amplitudes, while Ks of 
the order of magnitude of the sound speed, as expected 
for an acoustic wave, and Ki and K2 are of the order of 



the "spring constant" 



"(1.2) 



With these new notations we get: 



(I) 



9.(i?y') 



dAH'2^') 



dAH'J') 



Qsvl 



1 



Hl'^Hl'^ 



2S1/2 K1K2KS V^ 
uj„vl 1 



^1/2^1/2 



2S1/2 KlK2Ks ^Cg,Cg,Cg, 



El/2 K1K2KS ^c 
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r ig. 3. On this figure are presented two pairs of warps which can be 
coupled with the spiral. At each x, we see on the graph that CJ\+CJ2 — ^s^ 
and also that cD^/ + UJ21 — ujs- The grey band at the left arises from the 
dispersion relation of warps, which cannot propagate where \uj\ < /^. 



In order to be able to estimate the behavior of the 
solutions of system (I), we must analyze the geometry and 
localization of the coupling. 

4-. 2. Localization of the coupling 

To this point we have considered the coupling between 
three waves: the spiral wave, and two low amplitude warp 
waves which satisfy the selection rules. But any pair of 
771 = 1 warps obeying the frequency selection rule will 
also obey the other ones, leaving us potentially with a 
continuous set of warp pairs coupled to the spiral. The 
question is now to determine whether one or several cou- 
ples of warps can be preferentially amplified by the spiral 
wave, and if so, which one(s) and where. In order to clar- 
ify the situation, we present the problem graphically on 
figure H On this figure, each wave is presented by a curve 
(a segment in the shearing sheet approximation, but this 
does not lead to a loss of generality). For each curve, x 
is given by the implicit relation uj = uj — nin{x), where 
we have taken mg = 2 and toi — 1712 = 1. It is an easy 
matter to check that the slope of the curves describing the 
warps, at a given x, must be twice the slope of the curve 
describing the spiral. For the reason already explained in 
section 2.2.2 , if two warps satisfy the w-selection rule at 
some X, they satisfy it at any x. Since x has its origin at 
the corotation of the spiral, the curve of the spiral starts 
from (0,0). The warps cannot propagate in the vertical 
forbidden band [— /^,-|-/u], since their dispersion relation 
is (in the infinitely thin disk limit) w^ = ^"^ + 2nGYiq 
(taking into account the finite thickness of the disk and 
possible compressional effects maintains the existence of 
such a forbidden band, see e.g. Masset and Tagger 1995). 
For a given warp, the forbidden frequencies can be con- 
verted into a forbidden band in x as shown on figure for 
warp 2. 



We have depicted a narrow forbidden band since /x is 
expected to be small compared to k if the rotation curve 
is nearly flat. In the ideal case where the rotation curve is 
really flat, k — \/2^ giving ^i — Q. 

According to the expression of system (I), we see that 
the coupling can be strong where one or several group 
velocities vanish. Physically, this corresponds to the fact 
that the waves take a long time to propagate away from 
the region where coupling occurs, leaving them ample time 
to efficiently exchange energy and momentum. 

As already mentioned, we shall now use the group 
velocity of a stellar spiral for the spiral wave, although 
we keep our purely gaseous coupling coefficient. By do- 
ing so, we avoid the tedious handling of a mixed kinetic- 
hydrodynamic formalism without a great loss of precision. 
The group velocity of the spiral vanishes near the Lind- 
blad resonances (here at the Outer Lindblad Resonance, 
OLR) and the group velocity of the warp vanishes at the 
edge of the forbidden band, where an incident warp must 
be reflected Q 

For these reasons we have four qualitatively distinct 
cases to consider: 

Far from the OLR of the spiral: 

— We consider two warps, none of which close to its for- 
bidden band; 

— Or a pair of warps, one of which lies at the edge of its 
forbidden band. 

At the OLR of the spiral: 

— We consider two warps, none of which close to its for- 
bidden band; 

— Or a pair of warps, one of which lies at the edge of its 
forbidden band. 



4.2.1. Far from the OLR 

Let us analyze the first case. The graph is presented in 

figure 0. We choose a value of x such that neither warp 

1/2 
is close to its forbidden region. We define Ht by iJ/ = 

H^' + H2 . Summing the two first lines of system (I), 
and after some transformations, we obtain: 



dxHl'^ = 



^SJ^I 



2K1K2 y/c^^c^ S 



:H, 



1/2 



(4) 



Our purpose is to determine whether the coupling will 
be sufficient in these conditions to extract the warps from 
the noise level. As long as the warps remain at this low 
amplitude, their influence on the flux of the spiral is neg- 
ligible. We then assume that a/I] is constant. This natu- 
rally leads us to define an e-folding length for the warps. 

^ We can consider such a reflection as perfect, with equal 
incident and reflected fluxes, and with vanishingly small trans- 
mission of another warp wave beyond the forbidden band, see 
Masset and Tagger 1995. 
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h ig. 4. This figure presents the first case studied in the text. Wc 
analyze the efl^ciency of coupUng at the radial position given by the hor- 
izontal dotted line, far from the OLR of the spiral (top of the diagram). 
Both warps are far from their forbidden band. 



OLR 



Warp 1 



Warps 

Forbidden 

Band 




|i 



r ig. O. This figure shows the second case studied in the text. We 
analyze the efficiency of coupling at the radial position given by the hor- 
izontal dotted line, far from the OLR of the spiral (top of the diagram). 
Now one of the warps is located just at the edge of its forbidden band. 



If this length happens to be short (in a sense which will be 
defined below), the warps can be strongly amplified from 
the noise level. Here, the e-folding length is: 

■^e — z n 

iosv^ a 

Taking into account the expression of Ki given in ap- 
pendix y, and noting that the group velocity of warp i 
can be written as Cg. = an/QCoi^ we can rewrite this as: 




where we have taken the typical value (see the dia- 
gram) ujs — k/2. 

Now let us follow these two warps in their motion out- 
ward. After traveling over a distance Ae, their quadratic 
total flux Ht has been amplified by a factor e^ . Will ampli- 
fication at this rate be sufficient for the warps to extract 
a significant fraction of the fiux of the spiral ? For this to 
occur, the spiral and the beat-wave of the warps should 
be able to maintain a well-defined relative phase as they 
travel radially, lest the coupling term oscillates and gives 
alternatively a positive and negative energy flux from the 
spiral to the warps. We meet here a condition which is 
the WKB equivalent of the selection rule on the radial 
wavenumbers: since the WKB wavenumbers of the waves 
do not a priori obey this selection rule, they will decorre- 
late over a few radial wavelengths, k~^ r^ a/fl, where k^ is 
the radial wavenumber of any of the waves involved. Thus 
if kxXe is small the warps can be strongly amplified be- 
fore they decorrelate from the spiral, whereas if k^Xe is of 



the order of 1 or larger the coupling will not significantly 
affect them. Here, we have: 



f^x^e — 



4 Vi^ /wiWaV 



a 

~io 



Then the waves are in an adequate phase condition 
over one e-folding length at best, and we conclude that 
the warps cannot be extracted from the noise level in this 
first case. 

The second case, illustrated in figure ^, still corre- 
sponds to coupling far from the OLR of the spiral, but 
with one warp near the forbidden band. It is an easy mat- 
ter to check from the linear dispersion relation that the 
group velocity of a warp near its forbidden band is given 
by: 



Cg(Ax) 




(5) 



where Ax is the distance to the forbidden band. This 
formula has been obtained in the WKB limit, but we do 
not expect very different results in the general case. This 
leads us to replace in the expression of Ae a term of the 
type: Uj Jc^ by an integral of the type: 

^ d(Ax) 
3F(A^ 

We have adopted the length scale H because it corre- 
sponds to the range over which equation (m is valid, close 
to the forbidden band. 

The ratio of these two quantities is: 






1/4 



Q3/4 
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and is typically of the order of 1 (the exponent 1/4 
does not allow it to change very much). Furthermore, the 
coefficient ^Jujiuj^J^ is now about y^fi/Q,, i.e. about 
1/2. This is not sufficient to make the value of qX^ much 
smaller than 1. 

Thus we can conclude that far from its OLR the spiral 
cannot transfer its energy and angular momentum to the 
warps. 

4.2.2. Near the OLR 

We turn now the last two cases, corresponding to coupling 
close to the Outer Lindblad Resonance of the spiral. What 
does close mean here ? We said that we must study sep- 
arately the behavior of coupling near the OLR because 
the group velocity of the spiral in a stellar disk tends to 
zero as its approaches the resonance. In fact it does not 
really vanish (Mark, 1974), but we can in a first approx- 
imation consider that it does, with a decay length scale 
of the order of the disk thickness. Since this behavior is 
associated with the resonant absorption of the spiral wave 
at its OLR, by beating with the epicyclic motion of the 
stars, we expect that it should survive even in a realis- 
tic mixture of stars and gas. From the dispersion relation 
for a stellar spiral given by Toomre (1969), and using the 
asymptotic expression of Bessel functions, one can write 
the group velocity of the spiral near its OLR as: 



-gs 



"M^ 



3/2 / \ 3/2 



So that now the group velocity is a power 3/2 of the 
distance to OLR. 

We can slightly transform our expression of Ag, the 
e-folding length of the preceding section, as: 



A. 



4a 



V^ 



^^^"^Ks 



LU1UJ2- 



H 



,1/2 

1/2 "9S 



This expression is more appropriate in this new situa- 
tion, since (at least when the fluxes of the warps remain 

1 /2 

small compared to the flux of the spiral) Hg is constant 
(while S/cr is not) and now Cg^ varies. As before, we have 
to see how many e-folding lengths are contained in a scale 
length of the order of a/fl, where we can use the above 
expression for the spiral group velocity. Then we have to 
replace HjcJ^ in our previous estimate by L dx/cjs {x) , 
leading us to multiply the number of e-folding lengths over 
a scale H by: 



dx 



= 4 



Furthermore, we had taken in the previous section the 
mean value Cos = k/2. Now since we are close to the OLR 

we have to take ujs — k. Finally, with our typical figures, 

1/2 
instead of having only one e-folding length for H^ over 



a scale length H, we find eight such e-folding lengths, i.e. 
near the OLR the flux of the warps can be multiplied by 
g2x8 _ j^qT qj: course this flux cannot become arbitrarily 
large, so that when it reaches the spiral flux we can con- 
sider that most of the spiral flux has been transferred to 
the warps. Even though our factor of 10^ could be reduced 
by many effects, we consider that mode coupling near the 
OLR is efficient enough to allow the warps to grow from 
the noise level until they have absorbed a sizable fraction 
of the flux of energy and angular momentum carried by 
the spiral from the inner regions of the galactic disk. 

We now have to check in more details the effect of 
the localization of the coupling coefficient. Indeed from 
the system (I) we see that, if the amplitudes of the waves 
varied more rapidly than the coupling coefficient, their 
oscillations would reduce or even cancel over a given ra- 
dial interval the efficiency of the coupling. This must be 
considered because the same effect which gives a vanish- 
ing group velocity for the spiral near its OLR also causes 
their radial wavenumber to diverge. 

However, from the expression above for the group ve- 
locity of the spiral, we see that: 

so that 

K ^ [roLR -r)"^/^ . 

Now the phase of the spiral, in the WKB approximation, 
varies as 

exp (i / krdr) 

and the integral converges, so that in fact the phase tends 
to a flnite limit (instead of varying rapidly as one might 
think from the divergence oi kx)- Thus in the expressions 
on the RHS of system (I), the quantity which gives the 
dominant behaviour is indeed the (integrable) divergence 
of the coupling coefficient, in the vicinity of the OLR. This 
justifies that, in our estimates, we have integrated over the 
variations of this coefficient, keeping the other quantities 
approximately constant. It also justifies our statement, in 
section 2.2, that the generation of the warps is impulsive- 
like. 

This does not tell us yet which couple of warps is 
preferentially amplified at OLR, according to their fre- 
quencies. In fact the coupling coefficient still depends on 
these frequencies, through the term in \/ljiui2- This fac- 
tor, associated with the group velocities of the warps, fur- 
ther increases the efficiency for the most "external" pair: 
[ui = K—fJ., UJ2 = /i] (by a factor ~ 2 compared to the most 
central pair of warps if we assume /x/O ~ 1/4). This dif- 
ference is amplified geometrically at each e-folding length, 
and thus finally we can consider that the pair [/i, k — /i] 
is preferentially emitted at the OLR. This result is illus- 
trated in figure |6|. 
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r" ig. D. This figure shows the final result of our study. The coupling is 
strongly eflricicnt at the OLR, and only there, and results in a conversion 
of the spiral into two warp waves, at 0)2 — M and cDi — k — /^. 



Now we have determined how and where the spiral 
is "converted" into warps. But we still have to find how 
the energy is distributed between the warps, and in which 
direction they propagate. On figure |6| we have illustrated 
the warp 2 with an outward-pointing arrow, since we know 
that it is necessarily emitted outward due to the presence 
of its forbidden band. But warp 1 can consist of two waves, 
propagating inward as well as outward, and we must now 
consider the detailed balance of the energy transferred to 
the three warps. This can be done using the conservation 
of angular momentum. The relation between the energy 
flux H and angular momentum flux J is, for the warps as 
well as the spiral (CoUett and Lynden-Bell 1987): 



J 



mH 



since angular momentum and energy are advected at 
the same velocity Cg. Angular momentum conservation 
gives: 



(+) 



Js = jr + Ji 



r(-) 



J2 



where Js is the momentum flux lost by the spiral, while 
Ji, J2 are the fluxes gained by the warps, and the (+) et 
(— ) upper-scripts denote respectively the warps 1 emitted 
outward and inward. Writing Ji = J} + J} , and using 
the the fact that the waves 1 have the same frequency, 
we get Ji = niiHi/u} = Hi/uoi, giving from energy and 
momentum conservation the system of equations: 



Hi+H2^ Hs 

Hi/uJi + H2IUJ2 — 2Hs/uJs 



Defining /3i and (32 as the ratios Hi/Hs and H2/HS, 
we find (3i = LUi/uJs- 

Using the frequencies of the preferred pair of warps 
leads to: 



/3i = 

(32 = 



20 + K 

n + fi 



2n + K 



A strictly flat rotation curve would give: (3i = 0.71 
and (32 = 0.29, while a nearly flat rotation curve with 
^i/n = 0.25, would give: /3i = 0.63 and (32 = 0.37. 

Now it is an easy matter to check that the fluxes car- 
ried away by warps l'^' and 1*^"^ are equal. The easiest 
way to do it is to write for the pair of warps 1 a system 
similar to the one used above. This system is degenerate, 
but introducing a slight offset between the frequencies and 
making it tend to zero we find that they share equally the 
flux (3iHs. Physically, this corresponds to the fact that 
the source of the waves is "impulsional" (very localized), 
so that they are emitted equally without a preferred di- 
rection in space. 

From the numerical values given above, we can deduce 
that each of the three waves (two emitted, one reflected) 
carries away one third of the flux extracted from the spiral, 
as a very good approximation. 

4-3. Detailed physics at the OLR 

The situation exposed above is idealized. Two remarks 
modify the simple computations, but should not qualita- 
tively alter the main physical conclusions. 

— First, the group velocity of the spiral does not really 
vanish at the OLR. This has been investigated by Mark 
(1974), who showed that the group velocity of the spi- 
ral at the OLR was about a few kilometers per second 
(typically 10 times lower than the value it has far from 
the OLR). This leads us to replace the factor 4 in our 
computation by about ^ \/T0. Nevertheless, we still 
have a "large" number of e-folding lengths in the re- 
gion "just before" the OLR. 

— The spiral, as it approaches its OLR, is linearly 
damped through Landau effect. Mark (1974) gives an 
absorption length by Landau damping which has typ- 
ically the same order of magnitude as the e-folding 
length we find for non-linear coupling. The two pro- 
cesses can thus compete, the spiral losing part of it 
energy and momentum flux to the stars (heating them 
and forming an outer ring), and to non-linear ly gener- 
ated warps, at comparable rates. It is thus very likely 
that the dominant process (non-linear coupling or Lan- 
dau damping) will depend on the local characteristics 
of the galaxy under consideration. 
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F ig. T. Ttiis figure stiows the spiral wave and the three warp waves 
excited near its OLR by non-linear coupling. See the text for more de- 
tails. 



5. Comparison with observations 

In this section we derive the expected characteristics of the 
warps produced by the coupling process. Let us first sum- 
marize the main results obtained in the previous section. 
The reader can also refer to figure ^ where the situation 
is depicted. 

— The spiral wave travels outward until it reaches its 
OLR. Before it reaches it, the coupling with warps is 
too weak to affect the spiral. 

— When it reaches the OLR (more precisely within a 
distance of the order of a/V. from the OLR) the spi- 
ral slows down and is efficiently coupled with warps. 
This results in the "conversion" of the spiral into three 
warps which carry the incident flux of the spiral away 
from the coupling region: two warps at frequency n—fi, 
traveling respectively outward and inward, and the last 
warp at frequency /i, which is emitted outward. 

— The flux absorbed from the spiral is approximately 
equally shared by the three warps, for reasons of con- 
servation of angular momentum and energy. 

— In the following we will neglect the absorption of the 
spiral wave by Landau damping. Thus our results on 
the warps amplitudes will be optimistic, but still give 
results of the correct order of magnitude since the 
length scale of Landau damping and non-linear cou- 
pling are similar. 



flux extracted from the spiral (the rest of the flux of the 
spiral being absorbed by Landau damping, as discussed in 
section 4.3). 

In order of magnitude, the flux of the spiral is given 
by: 



Hs — CggEs ^ a x a Y, 



disk 



-'disk 



Similarly, the fluxes of the warps emitted outward are: 
Hw — Cg„Ew ^ a X uj^rT^outZ 

In these expressions, S^isfc is the unperturbed surface 
density of the stellar disk, where the spiral propagates, and 
Eout is the unperturbed surface density of the external HI 
disk, where the two emitted warps propagate. 

Now if we adopt a radial dependence of Y^out such as 
given by Shostak and Van der Kruit (1984), that we will 
very roughly describe by : 

Yout{r) = i;oexp[-(r- ro)/(i] 

where rg is the radius of the optical limit of the disk, d 
the decaying length of surface density (a sizable fraction 
of rg), we can deduce an approximate radial dependence 
of the deviation Z . We find : 



Z 
H 

-'disk 



-'disk 



Misk ^(r-ro)/d 



3Eo 



75 Mq.pc" 



and 

So ^ 10"^"part.cm~2 

(see e.g. Huchtmeier and Richter 1984 or Shostak and 
Van der Kruit 1984). 

and we take as a value for the amplitude of the spiral 
a /Yjdisk ~ 0.1 (but it could be larger than this widely 
adopted value, see Strom et al. 1976). This leads us to: 



Z/H^ le 



[{r-ro)/2d] 



5.1. Amplitude of the warp 

We will now derive a rough prediction of the warp dis- 
placement resulting from our mechanism. For each of the 
three warps we can write: 



-sw-^iv 



1 



gs 



Es 



where the factor 1/3 comes from the result (see section 
4.2.2) that the three warp waves share nearly equally the 



Of course T,out decreases progressively as the warps 
propagate outward, and Z accordingly grows to preserve 
the flux. 

The dependency found above is typical of what is ob- 
served in isolated spiral galaxies. It is noteworthy — and 
we believe that this has not yet been noticed — that the 
flux of the spiral waves is of the same order of magnitude 
as the flux of the warps. We consider this coincidence as a 
strong point in favor of our mechanism — independently 
of its details — , and a challenge to any alternate model. 
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5. 2. Expected characteristics of the corrugation 

The computations concerning the reflected warp are 
slightly different. One should remember that it propagates 
in a region where both stars and gas are present, a fact 
which cannot be neglected when discussing its character- 
istics. In a previous paper (Masset and Tagger, 1995) we 
discussed the propagation of warps in a two-fluid disk. We 
found two types of waves: the first type corresponds to 
waves where both fluids move essentially together, while 
the second type corresponds to waves where the fluids 
move relative to each other; assuming that the stellar disk 
is much warmer and much more massive than the gaseous 
one, in this type of wave the stars stay essentially motion- 
less while the gas moves in their potential well. We con- 
sider that this wave, which we call the corrugation mode, 
will be excited much more easily than the previous ones, 
since it moves much less mass and thus its amplitude for a 
given energy must be much higher. Its dispersion relation 
is: 



~2 2 

U = 11 



AirGp^ + 27rGSgg 



2 2 

:2«9'? 



where the subscripts g and * apply to gaseous and 
stellar values, and where p* is the stellar density at mid- 
plane. Because of the moderate thickness of the stellar 
disk, the square frequency v^ — AnGp^, is much larger 
than the other characteristic frequencies of the problem, 
and in particular larger than the frequency cj ~ k — /i 
of the wave which can be non-linearly excited. A solution 
can nevertheless be found because inside the OLR of the 
warp the last term in the dispersion relation is negative, 
and can for large enough values of q balance the dominant 
term i/f . Neglecting the self-gravity of the gas, and writ- 
ing (D^ <C 47rG'p* (a realistic assumption since the ratio of 
these quantities is about one tenth), we find: 



^/47rG'p» \/k? — uj^ 



Assuming a sech^{z/H) vertical stellar density profile 
{i.e. the disk is dominated by the local stellar gravity, 
which is a reasonable assumption), this can be rewritten 






(6) 



Taking a^ja^ ~ .2, and iJ* of the order of one kilopar- 
sec, we find an order of magnitude of one kiloparsec for the 
wavelength of the corrugation. This is consistent with the 
observed wavelength for the corrugation, in the Milky Way 
(Quiroga and Schlosser 1977, Spicker and Feitzinger 1986) 
or in other galaxies (Florido et al. 1991a). Despite the 
observations by Spicker and Feitzinger of a one kiloparsec 
corrugation in our galaxy, it seems that the typical wave- 
length of corrugation is rather about 2 kiloparsecs, or even 



larger (Florido et al. 1991b). The discrepancy of 2 be- 
tween observations and our estimate could be due either 
to the magnetic field which affects the gas motions and 
thus should imply the use of magnetosonic speeds instead 
of the acoustic speed, or to a lack of resolution of some 
observations. From equation (pi), we find that the wave- 
length of corrugations increases as we approach the edge 
of the optical disk, in good agreement with observations. 
Now let us give an order of magnitude of the expected 
amplitude of the corrugation. As already explained, we 
can write that its fiux is about one third of the spiral flux. 
Now the fiux of the corrugation is roughly given by (the 
index c is for corrugation): 



H,= 


^..E 


C 


where: 




duj 


O) K kH^, 


<'9. = 


dq 


~ q ~ q qH^ 


since gi?, ^ 1, and: 


E, = \vl^gZ^ 


Hence: 


^2 S.a^/ay 


wh 


ich gives, taking into a 


2, 





.•2/„2. 



Taking H^ ~ 1 kpc and cr/S* '^ .1 to .2, we obtain 
an expected amplitude for the corrugation of about one 
hundred parsecs, which is typically the observed value. 
Note that, as discussed above, energy can also be given to 
the other modes, involving both the gas and the stars, for 
the reflected warp. Their expected wavelength is larger 
than the disk radius, so that their line of nodes should 
be almost straight, corresponding to the observations of 
Briggs (1990), and so that the corrugation can easily be 
detected when superimposed on such large scale bending 
waves. 

5.3. Straight line of nodes 

We have found thus far that our mechanism leads to re- 
sults that fit well the observations for the amplitudes of the 
"outer" warps (the two waves emitted beyond the OLR of 
the spiral) and of the corrugation, and for the wavelength 
of the corrugation. We now turn to the wavelength of the 
outer warps. Actually the line of nodes of the observed 
warps is always nearly straight, leading to radial wave- 
lengths larger than the galactic radius. This is known as 
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the problem of the straight line of nodes, and is one of the 
most serious challenges to models of warps. In his work on 
the "Rules of Behavior of Galactic Warps", Briggs (1990) 
has emphasized this remark: the line of nodes of the out- 
side warps appears nearly straight or with a slight spiral- 
ity, always almost leading. Later observations have shown 
that outside warps could also be slightly trailing, but al- 
ways nearly straight. This can be interpreted from our 
mechanism, which predicts the emission of two warps out- 
ward from the OLR. In fact, the frequency of the first warp 
is fi, leading from the dispersion relation to a wavevector 
kx = 0, hence a straight line of node. The second one is 
emitted at a) ~ k — /i, i.e. near its OLR. From the disper- 
sion relation of warps taking into account the compress- 
ibility effects near the OLR (see e.g. Masset and Tagger 
1995), we see that k^ must tend to zero in order to bal- 
ance the behavior of the resonant denominator, thus also 
leading to a straight line of nodes. 

Of course this is a qualitative and simple interpreta- 
tion, since we use results obtained in the framework of the 
WKB assumption (tightly wound waves) where kx is pre- 
cisely not allowed to become small. Hence this discussion 
about the straight line of nodes is just the very beginning 
of a more complete work which should take into account 
non-WKB effects and the cylindrical geometry. This work 
is in progress, using the numerical code introduced in Mas- 
set and Tagger 1995. Finally, we should mention that the 
combination of two warps outside OLR can easily justify 
the large and arbitrary jump in the line of nodes across 
the Holmberg radius observed by Briggs, since the relative 
phase between warps 1 and 2 is arbitrary. 



6. Discussion 

6.1. Summary of the main results 

We have written the non-linear coupling coefficients be- 
tween a spiral and two warps, taking into account the 
finite thickness of the galactic disk and the three dimen- 
sional motion in a warp. The efficiency of this mechanism 
is too weak in the stellar disk, except at its outer edge 
where the spiral slows down near its OLR and can be effi- 
ciently coupled to warps and extract them from the noise 
level, transferring to them a sizable fraction of its energy 
and angular momentum (with the rest transferred to the 
stars by Landau damping). As a net result we can con- 
sider that the spiral is transformed into two warps. One 
of them receives a third of the energy of the spiral and 
can propagate only outward. The other receives the re- 
maining two thirds of the energy of the spiral, and splits 
into one wave traveling outward and one traveling inward, 
which we tentatively identify as the corrugation observed 
in many galaxies. 

Order of magnitude estimates of the respective ampli- 
tudes of the outer warp (in the HI layer) and of the inner 
one (the corrugation), as well as of their wavelengths, lead 



to values in good agreement with observations. In partic- 
ular we find from the observed amplitudes that the energy 
and angular momentum fiuxes of the spiral and the warps 
are comparable, a result which receives a natural expla- 
nation in our model. 

The expected long wavelength of the outer warp is in 
agreement with the observed trend to forming a straight 
line of nodes. Furthermore, our mechanism justifies the 
phase discontinuity of the warp at crossing the Holmberg 
radius (since the Holmberg radius is expected to be close 
to the Outer Lindblad Resonance of the spiral, where our 
mechanism takes place). 

6.2. What about the halo ? 

In the preceding sections we have not mentioned the role 
of a massive halo. We will not discuss a misalignment be- 
tween the disk and the principal plane of the halo, since 
such a misalignment was introduced in an ad hoc model 
by Sparke and Casertano (1984) in order to explain the 
continuous excitation of the warp in isolated galaxies. Our 
mechanism does not require such an hypothesis. 

A typical halo is about ten times more massive than 
the disk. If this halo has a spherical symmetry, it does not 
affect the vertical restoring force and thus leaves our re- 
sults unchanged. On the other hand a flattened halo would 
give rise to a strong restoring force for a test particle (or 
a disk) displaced from the equatorial plane. The vertical 
oscillation frequency /i would be about three times the 
epicyclic frequency. Hence the forbidden band of figures 
H-o would invade more than the whole disk. Our analy- 
sis would then have to be altered, but we believe that its 
main result would remain valid. 

Indeed, in order to find two bending waves that would 
propagate in the vicinity of the OLR of the spiral, we 
might invoke the "slow" branch of the warp dispersion re- 
lation, which uses intensively the compressional behavior 
of the gas (let us emphasize that we have already made 
use of this slow branch to find the wavelength of corruga- 
tion, which is in the same manner a bending wave of the 
gas in the deep potential well of the stars). This would not 
strongly change the expected efficiency of our mechanism, 
and would even rather increase it slightly, since we would 
be dealing with slower modes. Another solution would be 
to consider the pair of warps [— /z, k + n], which are both 
simultaneously propagating (with the former being retro- 
grade) . The efficiency of coupling with such a pair would 
be, in our framework (shearing sheet and steady state) ex- 
actly the same as the one we derived for the pair [/i, k — /i] . 
Now the details of the scenario would change a little. We 
would have only one wave emitted outward and two in- 
wards, but the order of magnitude of their amplitudes 
would not change. Note that observations of corrugation 
in external galaxies (Florido et al. 1991a) seem to reveal 
two distinct wavelength for corrugation waves. 
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Another matter concerns the response of the halo to 
the warp. A fully consistent treatment of this problem 
should be numerical. However, it is possible to get ana- 
lytical trends by considering the halo as a thicker (and 
then warmer) disk than the stellar/gaseous disk. This can 
be done using the two-fluid dispersion relation of warps 
(see Masset and Tagger 1995). The two-fluid dispersion 
relation reveals three modes: 

— A mode which involves essentially the thinner and 
lighter component, leaving the other one (the halo) 
nearly motionless. This mode is the one described in 
the preceding lines, and corresponds to the gaseous 
corrugation within the stellar disk as well as to the 
compressional mode of the whole disk (stars plus gas) 
in the rest potential of the halo. 

— The classical two modes (slow and fast) of the one- 
fluid analysis, which involve similar motions of both 
fluids. This means that the halo would participate in 
the different warps, as well as in the spiral, and for- 
mally the problem would be the same as the one de- 
scribed throughout this paper, with the difference that 
we would have to use the total surface density, involv- 
ing the projected one of the halo, and the total thick- 
ness, which is the halo thickness, of the order of r/2. 
Each of these quantities should then be multiplied by 
about ten. Hence the test particle frequency vz would 
remain the same, since it is proportional to yJH/H . We 
then see on equation (H) that each factor would remain 
unchanged, except the group velocities whose expres- 
sion would probably be slightly different, but could not 
differ much from the sound speed; hence the e-folding 
length would have the same behavior and order of mag- 
nitude. Furthermore, a participating halo would give a 
more efficient and simple explanation to the problem 
of the straight line of nodes, since the wavelength of 
warps, for a given excitation frequency, increases with 
the surface density. 

Figure || shows that the participation of the halo to 
the global modes (i.e. the two waves corresponding to 
the one-fluid ones) is as strong as or stronger than that 
of the disk. Furthermore the wavelengths of these modes 
agree with the observations of a straight line of nodes. 
However, this figure has been obtained in the WKB as- 
sumption and moderately thick disk formalism (Masset 
and Tagger 1995), which are both far too rough in this 
case, so that a numerical solution taking into account the 
halo response would be necessary. 

In the same manner, if we want to estimate the cou- 
pling efficiency taking into account the halo, a problem is a 
determination of a /H for the spiral involving a participat- 
ing halo. If we assume a behavior of matter independent of 
z, then a/Yj should still be about 0.2, but once again only 
a numerical solution would allow us to give the real verti- 
cal profile of perturbed quantities in the spiral and in the 
warps. Furthermore, the problem is much more complex 




r ig. 8. In this figure we present the three modes for bending waves in 
a two-fluids system (disk -denoted by a star— and halo —denoted by a h— 
). We have adopted the following parameters: Yih/Yi^ — 10, ^ / k — 0.1, 
lj/k — 0.75. The results (the elongation Zf^ of the halo over the elon- 
gation Z^ of the disk, and the wavelength) are presented as a function 
of the ratio a^/af^. For a typical galaxy, this ratio should be about 0.2. 
The solid line represents the classical infinitely thin onc-fiuid mode. We 
see that its wavelength is far greater than the disk radius (hence it could 
be a beginning of the justification of the straight line of nodes prob- 
lem) and that the disk and halo have nearly the same elongation. The 
dotted line represents the onc-fiuid compressional mode. Its wavelength 
is about four times the disk radius, and we see that the halo is more 
warped than the disk since the halo elongation is about 1.6 that of the 
disk. Finally the dashed curve represents the corrugation-like wave, i.e. 
here the "corrugation" of the whole disk (stars and gas) in the rest po- 
tential of the halo, which appears to be in that case nearly motionless ( 
Zh « Z,). 



if we want to treat it fully correctly, since for distances 
to the galactic center less than the Holmberg radius we 
have to cope with three fluids: the halo, the stars and the 
gas, so the method which consists in treating the disk as 
a unique fluid will give only approximate results. 

Nevertheless, despite the large number of allowed 
modes of propagation of bending waves in this realistic 
case, we see that global modes (i.e. involving a motion of 
the halo) will lead to the same efflciency as found in this 
paper. Numerical simulations would be needed to give an 
accurate evaluation of the c-folding length and to deter- 
mine the fraction of incident energy that each wave carries 
away. 

Finally, let us mention that the halo should be a col- 
lisionless fluid {i.e. star-like and not gas- like) if we still 
want the spiral to slow down at the OLR. Observations 
suggesting that the halo or outer bulge is composed of 
very low mass stars would conflrm this hypothesis (see 
e.g. Lequeux et al. 1995). 
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6.3. Suggested observations for checking the spiral- warps 
coupling mechanism 

In this paper we have shown that the spiral wave is con- 
verted into bending waves observed as an outer HI warp 
and a corrugation wave. We have aheady emphasized the 
observational fact that the fluxes of these three waves have 
the same order of magnitude for a typical galaxy. It would 
be interesting to test this coincidence with a better accu- 
racy, and to search for a correlation between the fluxes of 
the spiral, the warp and the corrugation. The fluxes can be 
straightforwardly obtained from such observables as the 
disk surface density, the rotation curve (which gives val- 
ues for K and /x), the disk thickness, the perturbed surface 
density of the spiral (see e.g. Strom et al. 1976), and the 
deviation of the disk to the equatorial plane — although 
some of these quantities are better measured in face-on 
galaxies, and others in edge-on ones, making a complete 
set of observations in an individual galaxy a challenge. A 
correlation between the fluxes for a significant number of 
galaxies would be a strong argument for our mechanism. 

Appendix A: Formal derivation of the coupling 
coefRcient 

A.l. Derivation without the shear terms 

We use the hydrodynamical equations (continuity and Eu- 
ler) and the Poisson equation. We write the first two equa- 
tions in tensor formalism, in order to avoid too lengthy 
expressions. The continuity equation is: 



dtp + d'{pVi)^0 



where the vector Vi (i = 1, 2, 3) is a tensor which rep- 
resents (f7, y, VF), and where we implicitly sum over re- 
peated indices (Einstein's convention). 

The Euler equation is: 



dt 



, dt[po{l + s)] 
Po(l + s) 



'9V(A.2) 



In this expression the last term of the L.H.S. represents 
the Coriolis acceleration, and e^^'^ the Levi-Civita symbol. 
We also introduce a convenient notation, s = p/po- The 
reader can note that, in the convective derivative {V.W)V, 
we have not written the shear terms. We neglect them in 
this first approach, and we will check their influence at the 
end of this section. 

In the Euler equation the pressure term can be ex- 
panded to second order as: 



(1 - s)d\s + 



Po 



Finally, the Poisson equation is: 
Act) = 47rGp 



We want to know the energy injection rate into warp 
1 from the spiral and the warp 2. In the linear regime, we 
would have: dtEi = 0. Here, we will have a non-vanishing 
expression due to the presence of the spiral and warp 2. 

Since we treat the warp 1, we expand the hydrody- 
namic equations to second order in the perturbed quanti- 
ties and project them onto e~*'^^*. 



For the continuity equation (A.l) this gives 



+ d,{poVI) + d,{psVn + d,ip;V^) = 



dpi 
dt 



and for the Euler equation ( |A.2| ): 
-a^d'si - 9Vi + a'^ssd's* + a^s^d'ss 



(A.4) 



(A.5) 



Our purpose is to write the temporal derivative of an 
expression which can be interpreted as the energy density 
of the warp (energy density per unit surface) . For this we 



multiply equation (A.5) by c?p\l pa and equation (A.4) 



by po: we integrate over 
conjugates. This gives: 

dt [[po{V^V*+a^sisl)] + 



add them and their complex 



po"l^iii9Vi + c.c. 



(A.6) 



a^poV*d'si- / a'sid^ipoV^)- / paVlV^djV, 



(A.l) - / PoViyfdjV^ + J Poa'V^^ssd's; + J poa^V^s^d'ss 



a^sld^ipsV^ 



a^sld^{plV^) + c.c. 



In this tensorial notation we will freely use the rule of 
integration by part {i.e. transpose derivative operators 
with a change of sign) in the linear terms of this equation: 
indeed if the index i is z, this is really an integration by 
parts; on the other hand Hi = x or y the derivative is just 
a multiplication by iki (in our WKB analysis), so that e.g. 



sid,{poVt 



-poVl*d^si 



(A.3) 



Furthermore, we must remember that the integration is 
only over z, so that we do not need to worry about bound- 
ary terms the integration by parts would introduce in an 
integration over x, and which would represent energy and 
momentum flux at the radial boundaries of the integration 
range. 

However, we cannot blindly apply this simple formal 
integration by parts on the non-linear terms. A careful 
look at the details of the integration by parts for each value 
of the index i is necessary. For the dimensions {i = y and 
i = z) the formal integration by parts would give correct 
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results, even on the non- linear terms. But for i = x, an 
integration by part would imply a selection rule of the 
type kxs — kxi + kx2, and we have intentionally avoided 
such an hypothesis. Since we want to stay, for simplicity, in 
the framework of tensorial formalism, we have to perform 
transformations that are correct whatever the index i and 
we do not integrate by parts the non-linear terms in this 
first part of the derivation. 

Equation (A6) already has the form we want : the first 
term in its L.H.S. is the temporal derivative of the kinetic 
and internal energy of warp 1. We can note that the term 
linked to Coriolis acceleration has vanished, since the Cori- 
olis force is always normal to the velocity and thus does 
not work. 

The second term in the L.H.S. must be rewritten to 
be interpreted as the temporal derivative of the potential 
energ y. We integrate it by parts and then use equation 
(A. 4) to transform the result. We can also note that the 
first two terms of the R.H.S vanish together (integrating 
one of them by parts) . We obtain: 



(whether the integral extends only over z or over the 
whole space, since the integrated quantity does not depend 
on X or y). 

Now: 

5* dtpw = dt{S ^ pw) 

since S is time independent. This is linked to the use 
of the Poisson equation, i.e. to the fact that the potential 
propagates instantaneously. Then one can write: 



/ (l)idtp*i = / p*idr 



Thus we see that we can write this term as a potential 
energy associated to the warp, and this gives: 



dt / Po 



{VIV^^ + a'sisl) + ( is*0i + c.c. 



(A.8) 



dt [po{VlV*+a^sisl)] 



4>idtPi +C.C. 



(A.7) 



PQV*yp,vr- / PoV{yfd,v^+ / poa'v,\.ssd 



+ J poa^V,\s;d'ss - J a^sld,ipsVr) - J a^sld,ip;V^) 

+ J <l>*A{psvn + J ^*Aip;v^) + c.c. 

The second term in the L.H.S still does not appear as 
the time derivative of a potential energy. We have to use 
the Poisson equation, which we have not used yet. 

If we write S{x'^ —Xq) for the potential created at point 
a;* by a unit mass located at point Xq, then: 

(j) — S * p 

where * represents the convolution operator. 
Thus, 4)1 = S * pi, 0s = S" * P5, etc. . . . 
We write £{p,p') the form: 



J poVl VidjVr - J PoVl VfdjV^ + J poaVi* ssd'i 



+ / poa'Vls*^d^ss~ / a'sld,{psVr)- / a'sld^ip^V^) 



rAiPsv^ 



r^^^ip;v^) 



c{p,p') 



p*{r)[S*p']{r)cPr 



where the index E indicates that the integral extends 
over the whole space. Using the fact that S* is a real func- 
tion with spherical symmetry, it is an easy matter to check 
that (by writing explicitly the convolution and exchanging 
the integrals): 

c{p,p')^c{p',pr 

Hence £ is a variational form, and one can write the 
following equalities: 



K^tPi = C{dtpi,pi) ^ C*{pi,dtpi) 



Hence we have expressed the temporal derivative of the 
energy density (kinetic plus internal plus gravitational) Ei 
of warp 1 as a sum of integrals which involve the ampli- 
tude of each warp and of the spiral. In particular, one 
may note that in the linear case the energy of the warp 
is conserved. One important remark concerning the linear 
case is that, if we neglect the coupling term, we should not 
have dtEi — but dtEi + dx{cgEi) ~ since energy is 
transported at the group velocity Cg. Here the convective 
term has disappeared in the integration by parts of one of 
the terms in the RHS of equation (A6): for simplicity we 
have neglected the integrated terms, assuming periodicity 
in X. From a more complete derivation we would recover 
a linear term dx{CgEi), without affecting the non-linear 
ones. We will thus re-write (A8) as: 

dtEi + dxicgEi) = —El = Coupling Term 

A. 2. Effect of the shear term 

As emphasized above, we did not take into account the 
shear term when writing down the Euler equation. Writing 
this term can a priori modify the warp energy and the 
coupling term. 

Formally, in tensorial notation, the shear can be writ- 
ten as: 

V^ = 2A'^xt 
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where the matrix A is, in the frame xyz: 




With these notations, the convective term of the Eulcr 
equation reads: 

i2Ai''xk + V^)dj{2A'^xi + V) 

which can be rewritten as: 



^A''A)xk 



2A)V' 



2A^''xkdjV' + V^djV' 



The first of these four terms, proportional to A , is null 
(this is generally the case when the velocity is constant 
along a current line). 

The last term is the one we had without shear. The in- 
fluence of this term on the coupling coefficient has already 
been derived in the previous section. 

The intermediate terms are linear, and thus have no 
incidence on the coupling term. 

Hence, the temporal derivative of the energy density 
of the warp is modified by the addition of the two cor- 
responding terms, obtained by multiplying by V* and 
adding the complex conjugates. This gives: 



dtE + m[ / paV*A)v^] + m[ / poV*A^''xkd,v 



dtE 



where 9? stands for the real part. For the waves under 
consideration both of these terms are null: 

— The first one because, in the WKB approximation, the 
perturbed horizontal velocities U and V are in quadra- 
ture (epicyclic motion). Thus this term, given the ex- 
pression of A, is proportional to UV* , and therefore 
purely imaginary. 

— The second one because it is related to the azimuthal 
derivative of terms related to the warp energy density, 
and thus vanishes. 

Hence the addition of shear does not modify the cou- 
pling equation. 

Appendix B: Analytical computation of the 
coupling term 

B.l. Approximation of the eigenf unctions 

In the framework of our expansion to second order in the 
perturbed quantities (assuming weak non-linearities), the 
coupling integrals derived in the preceding appendix in- 
volve the linear eigenfunctions associated with the spi- 
ral and the warp. These eigenfunctions consist in the am- 
plitude and phase relation between the perturbed quan- 
tities (velocities, density, potential...), and their spatial 
variations. An exact, analytical or numerical, knowledge 



of these eigenfunctions is beyond the scope of this work, 
and in fact would add little to it since we are not inter- 
ested here in deriving detailed numbers appropriate to a 
given galaxy model but rather in the physics of non-linear 
mode coupling. We will thus use approximate expressions, 
which allow an easier access to this physics. Furthermore, 
our results are given in terms of variational forms. These 
are well known to preserve the important invariants in the 
problem (the energy and action densities) , and to be good 
estimates even with poor approximations to the eigenvec- 
tors. 

Thus we remain in the WKB formalism, and make the 
following assumptions: 

— We assume that the motion in the spiral is purely hor- 
izontal and independent of z. 

— We assume that the vertical velocity in warps is inde- 
pendent of z, and that the perturbed density is: 

Pw = —ZwdzPo 

From these assumptions we can deduce the perturbed 
quantities relative to a warp. We focus hereafter on the 
warp 2, but the results would evidently apply also to warp 
1. 

One gets U2 et V2, the horizontal components of the 
perturbed velocity, from the relation: 



kxM2 + ky^V2 = - 



q2^2 



62 + a^S2) 



Using the WKB hypothesis {ky^ <^ k^^ — 92), one 
obtains: 



t/2 = - 



92^2 



+ a^S2) 



V2^ 



2in 



q2LU2 



+ a^S2) 



where 4)2 is known from p2- 

It is noteworthy that these quantities are only a rough 
approximation to the eigenfunction of a warp; in partic- 
ular they do not fulfill the continuity equation. On the 
other hand our reason for taking them into account is that 
near the Lindblad resonances, where horizontal motions 
are large, we suspect that they might strongly contribute 
to the coupling terms. We will find later that this is not 
the case, so that forgetting these terms altogether would 
not change the result to leading order. 

The perturbed potential of the warp is given by the 
expansion to lowest order in qH of the expression given 
by the Green functions (see Masset and Tagger 1995): 

One obtains : 



Mz) 



-47rGZ2 



po{z')dz' 
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and then the other perturbed quantities associated We will make intensive use of integrals of the type: 

with the warp can be expressed, as functions of Z2 ex- 

clusively: /^ = / "^^i^^^dz 

J -00 Po 

U2 = —rr^ tt\4:TtG / po(z')dz' + a'^— — ] Z2 and 

^2-'^2\ Jo Pa J 

2in q2U2 ( , ^ r , ,., , ^ 2dzPa\r, ^2 = 

1/2 = -To o 47rG / po(z )dz + a Z2 

« ^2 - '«2 V h Po / 

All the other integrals are either straightforward or can 

W2 — ~iuo2Z2 be deduced from Ji or J2. 

For the evaluation of Ii , we use the hydrostatic equi- 

3„ — ^2 ^ librium in the unperturbed state: 



dzPo I Po{z')dz' 
Jo 



dz 



Po 



a dzpo = ~podz(t>o 



In the same manner, with the hypotheses mentioned 
above, we obtain for the spiral: Hence 

ujs a /■+°° 

Vs^ — :^ /i = - / d^(t)od;,podz 

IS ^ J -00 

2iVl LOg a Integrating by parts, one obtains: 



2tiG(j 

<Ps = -- 



/i = / Pod] 



3 



Qs 



Now 

a^ 00 = A0O- Ar0O 



E where A^^o, the radial laplacian of 00 is equal to — /i^ 

(Hunter and Toomre, 1969). 
Rewritmg the coupling term obtained in the previous ^^.^^ ^^^ p^-^^^^ equation, A0o = ^nGpo, one ob- 



section, and expanding the sums over repeated indices, we 
get: 



tams: 

h ^ 4:nGY.2 + ^l^T. 

where S2 is defined as: 
i{ls - q2)poUsS2<t'i + i{q2 - qs)poa^s*s{UiS2 - U2S1) /•+°o 



—El = / i{qs ~ q2)po(t)iU2S*s + PQSgW2dz(j)i (B.l) 



+PQa'^s*s{WidzS2 + W2dzSi) + iiqs - q2)poa^siS2U*s 

+ipoUiU2U*{qs - q2) - Poiq2ViV2U*s - Poiq2WiW2U*s 

+Pf)iqsU2ViVg + c.c. 

The derivation needs a frequent use of the following 
The final step consists in replacing the quantities in- integral: 
volved in these integrals by the approximate perturbed 



Po{z) dz 

Using an integration by parts, one can easily derive: 
I2 ~ — ^2 



quantities derived above. 

B.2. Computation of the coupling integrals 

In this appendix we compute the various terms in the RHS 



h^ I Po 



AtiG / pQ{z')dz' + a 



2^zP0 

Po 

Expanding the square in the integral: 

2 



2 
dz 



r— 1, f-rOO / fZ \ ^ 

of the coupling equation (^). j^ ^ IGir^G^pol po{z')dz'] +a^h+87rGa^l2 

We will not write down all the computations, which J-00 \Jo ) 

are lengthy and tedious. Let us just notice that the spiral 

eigenfunctions do not depend on z, so that each integral '"^'^ deauce. 

over z involves in fact only two eigenfunctions of warps 1 ^_2/^2y'3 

and 2. ^3 = 7, ^i^Ga^Y.2 + P^a^S 
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In order to express the coupling term , we group the 
terms as follow: 



d 
dt 



Ei= ijqs - q2)pQ(f>iU2S*s^ + ipoSgW2d^(j>i^ (B.2) 

+ Kqs - q2)PoUgS24>i +i{qs - q2)pQa^s*s{UiS2 - U2S1) 

+ ipoa^Sg{WidzS2 + W2dzSi)+i{qs - q2)poa^siS2Ul 
-■ , ' ^^ , ' 

+ ipoUiU2U*{qs - (72) -P0'iq2ViV2Us - poiq2WiW2U*s 

^ ^ '^ V '^ V ' 

Tr Ta T9 

+ paiqsU2ViVs +c.c. 

and we find for these terms the expressions: 

{qs - 92)92^2 



Ti 



CSi, — K^ 



4^^G^I]^ . ^ 2^2 



Z\Z20 



T2 = 47rGa;2 — Z\Z20 

rr . r^~ Is - q2^2 ry ry 



qs 



Ti = a'p.'{qs^q2) 



qiuji 



92^2 



Z\Z20 



T. 



n = 



AttG— + p^ ) (UJI + Cj2)ZiZ2(T 



{qs - 92)^5 



qs 



^ qs-q2 ~ 

T7 = UJs—2 



47tG— + p^j: ]ZiZ2(j 

qiq2UJiUJ2 



qs 



{Col ~ '^^)('^2 ~ '*^) 



Tg = — ojs 



'in'G^j:^ 2^2 22, 

X I AnGa^^ + p^a^ ] ZiZ2a 

52 - 4ri^ qiq2U}iU}2 



qs "^ K^ ((Df — K2)(a)| — K^) 



(Aii'^G'^Y? ^_^ 2^2 22 






- AnGa^-^ + ^^a^ ) ZiZ2(t 



n 



Ll!SLl!iLlI2ZiZ2(J 

qs 



4^2 



qiq2UJiLU2 



K^ {Ql — K^){U}2 — ^2) 



/ 4^2g2^2 

V 3 



- AnGa^— + p^a^ ] ZiZ2(t 

2-i 



waves were such that this product is real the coupling term 
(once added to its complex conjugate) would exactly van- 
ish, while it would be maximized by waves in quadrature. 
This variation with the phases is associated with the com- 
plex non-linear behaviors mentioned earlier, and we will 
not discuss it here (technically the phases can evolve non- 
linearly on the "slow" time scale of the mode growth and 
non-linear evolution, as compared with the "fast" time 
scale of the linear frequency). Hereafter we will for sim- 
plicity assume that the non-linear process has picked up 
waves from the background noise, or made them evolve, 
such that their phases maximize the coupling, i.e. Z\Z20 
imaginary. Thus we will from now only consider real quan- 
tities in the coupling equation. 

Factorizing X]i=i ^i, we obtain the following coupling 
equation: 



dt 



E, ^ 



peie2 — ^CJ5 



4^2 

k2 



1- 



417^ 



q2_ 
qs 



(B.3) 



+ (a;i(D2 +r] ~ 2fi ) — ljs + {p - p )[ UJ2 + —l^s 

qs \ qs 



+ ^^^[g2e2(/3 + a2^2) _ q^e^„2^2] 
K 



where we have introduced: 



Z\Z20 



We note that the coefficients in eq. (B2) are all imag- 
inary, and that all the expressions above give a real fac- 
tor times ZxZ2(y. Thus if the relative phases of the three 



r/^ = p'- - 47rGI]2/S 

13 = 47r2G2s2/3 + 02772 

/+00 
PQ{zfdz 
-00 

Q\ = a)iK/(a)]^ — K ) 
62 = a;2K/('^2 - >^) 



We clearly see that we are concerned with inhomoge- 
neous coupling since the coefficient depends on x, which 
implicitly comes from the w. In particular, the dimension- 
less coefficients 9i and 82 are of the order of unity, except 
near the Lindblad resonances of the corresponding warps 
where they can take very large values. 

Appendix C: Simplification of the coupling term 

In this appendix we simplify the expressions obtained in 
appendix to get an estimate of the coupling coefficient, 
i.e. of the efficiency of non- linear coupling. First we esti- 
mate the constants /3, 772 , and S2 which appear in equation 
(B.3). /3 can be written as /3/S (see appendix Q). From 
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the hydrostatic equilibrium of the unperturbed state we 
find: 



First, we compare v^^ and /i^. We have: 



PO JO 



47rG / pQ{z')dz' + c? 



AnGpo — d1 <j)Q dz' — — /i^z 






This gives from the definition of /s : 

/ + 00 
poM z dz ~ n YiH 
-CXD 

and thus: 



/i/2 



Now we have to compute 47rGS]2/S. One can say that 
E2/E is of the order of magnitude of pm, which is the 
density in the mid-plane of the disk. Now AnCpm is the 
frequency of vertical oscillations of a test particle in the 
rest potential of the galactic disk. We note this frequency 
Vz, where the z index reminds that this frequency depends 
on the vertical excursion of the test particle. Hence: 
222 

Finally we consider 8 = kCj/{ijP — k^). 8 is of the or- 
der of 1, except in the vicinity of the Lindblad resonances, 
where it become large. This is essentially the same ef- 
fect that was found in previous works (Tagger et al. , 
1987; Sygnet et al. , 1988) to make non-linear coupling 
very efficient near the Lindblad resonances of the coupled 
waves (although in that case the coupling coefficient was 
obtained in a kinetic rather than fluid description, i.e. 
adapted to the stellar rather than gaseous component of 
the disk). On the other hand, our assumption of weak non- 
linearities certainly breaks down if 8 becomes too large, 
i.e. exactly at the Lindblad resonances. A convenient esti- 
mate is that this break-down occurs at one sound-crossing 
time a/ K (or equivalently one epicyclic radius) from the 
resonances, giving an upper limit for 8: 

«2 1 



8 



max o 9 



g2ij2 



10 to 100 



Now we can estimate the different terms of equation 



B.3. The first term can be estimated as follows: 



— The term in inner brackets is of the order of 2. This 
estimate is of course expected to vary of a factor 2 or 3, 
but no more. In particular, it is irrelevant to speculate 
on a possible resonant effect due to a vanishing value 
of qs'. as we will find below, the spiral is tightly wound 
in the region of strong coupling. 

— It is an easy matter to see that 8182 is always about 
0.5, even in the vicinity of the Lindblad resonances, 
since d)i + 0)2 = ujs- 

— Finally, j3 is of the order of fi'^H^, as seen above. 

One can then deduce the order of magnitude of the 
first term 71: 



Ti = /iJ 



2 9192- 



-U!S 



(C.l) 



Let us evaluate now the second term T2 : 



4:TrGp„ 
GM/r 



4r 
ll 



40 



Now 51 is larger than p in the region of the galaxy 
where the rotation curve is nearly flat, hence ff 3> /i^. 
On the other hand, since uj r^ fl, we have v^ 3> 0)10)2 ■ 
Finally, we assume that (72 < Is (this is a good ap- 
proximation in the region of strong coupling; we will 
emphasize it again when we study the localization of 
the coupling ). 



We deduce that: 



T2 = V^LJS 



(C.2) 



Let us note that T1/T2, from eqs. (|CJj) and ( |C.2| ), is 
of the order of: 

where each factor is smaller or much smaller than 
unity. Hence the flrst term is always negligible compared 
to the second one. 

Let us now find an estimate for Tjj . We will derive it in 
the vicinity of a Lindblad resonance, maximizing 8. We 
note that: 

13 P^ , 

so that /3 -I- a^jU^ ~ a'^fj,'^. 
We estimate: 



% 



1 



1 



'3 ^ Qs 22 
— ~ —a p 92^772^"^ 
T2 K q^Hi vjus 




Hence TJj is smaller or comparable with T2 in the vicin- 
ity of a Lindblad resonance. Away from the resonances, 
the ratio 73/72 is still lower. Thus we deduce, since 7i is 
always negligible compared to 72, that the coupling term 
is always of the order of 72 : 



dt 



El ^ v^LbsZiZ2(y 



Appendix D: Expression of energies 

The coupling equation written in the previous appendix 
involves both the energy and the amplitude of warp 1. 
They are actually linked by an expression we wish to de- 
rive. 

The energy Ei reads (c/. eq. (|l|)): 



-El = 



po(^)(|(7iP + \Vi\'' + \Wi\^ + a^lsip + ^sl<Pi])dz 
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^z and 




r" ig. 13.1. This figure shows the velocity field of a 
warp aw^ay from the Lindblad resonance. 




1" ig. 13.2. This figure shows the velocity field, 
which becomes nearly horizontal, of a warp near a 
Lindblad resonance. 



After a straightforward calculation, in particular mak- 
ing use of the integrals of appendix H, one obtains: 



El 



1 



^2" 






~ 2 

■UJ^ 



M" 



SZ^ 



Far from the Lindblad resonances, it is easy to see that 
the first term of this expression of energy is negligible com- 
pared to the third, and a fortiori to the second. However, 
at the Lindblad resonances of the warp, the first term can 
become important and dominate the others. The physi- 
cal interpretation is that near the Lindblad resonance the 
kinetic energy associated with horizontal motions, due to 
the compressibility of the gas, becomes dominant. 

Hence, at the Lindblad resonance, on can have energy 
"hidden" in the horizontal motions associated with the 
warp, i.e. a large energy with a small vertical displace- 
ment. This is illustrated in figures D.l and D.2. 



On the other hand, for a spiral, the vertical motion 
never dominates even when compressibility becomes im- 
portant. This can be directly seen from the expression of 
the energy of the spiral: 



Es 



I- 



T2~ 



9| 






which does not show any resonant term. Thus the per- 
turbed surface density a represents fairly well the energy 
of the spiral wave. There are no hidden motions (hidden 
in the sense that they don't have incidence on the observ- 
able) similar to the hidden horizontal motions in warps. 
Thus we have: 

1/2 

^+^ L n ^ n I I, n n 



K,= 



4f^' 2„2o2M* 



~2 , 2 



for i=l or 2 



Ks^ 



^s 
ll 



4rj^ 
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